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ABSTRACT 

We compile black hole (BH) masses for ^ 60,000 quasars in the redshift range 0.1 < z < 4.5 included in 
the Fifth Data Release of the Sloan Digital Sky Survey (SDSS), using virial BH mass estimators based on 
the H/3, Mgii, and Civ emission lines. We find that: (1) within our sample, the widths of the three lines 
follow log-normal distributions, with means and dispersions that do not depend strongly on luminosity or 
redshift; (2) the Mgii- and H/3-estimated BH masses are consistent with one another; and (3) the Civ BH mass 
estimator may be more severely affected by a disk wind component than the Mgii and H/3 estimators, giving 
a positive bias in mass correlated with the Civ-Mgii blueshift. Most SDSS quasars have virial BH masses in 
the range 10^ - 10'" Mq. There is a clear upper mass limit of ^ 10'" Mq for active BHs at z > 2, decreasing at 
lower redshifts. Making the reasonable assumptions that the underlying BH mass distribution decreases with 
mass and that the Eddington ratio distribution at fixed true BH mass has non-zero width, we show that the 
measured virial BH mass distribution and Eddington ratio distribution within finite luminosity bins are subject 
to Malmquist bias if the scatter in luminosity at fixed true BH mass is uncorrected with the scatter in line width. 
Given the current versions of virial calibrations and their uncertainties, we present a model which reproduces 
the observed virial mass distribution, quasar luminosity function, and line width distribution of our sample; it 
has an underlying BH mass distribution which is a power-law with slope ^ -2.6, and a true Eddington ratio 
distribution at fixed BH mass which is a log-normal with mean dependent on BH mass 10"''^ for typical 
10^ Mq BHs) and with dispersion 0.4 dex. In this model, the observed virial mass distribution for the SDSS 
sample is biased high by ^ 0.6 dex within finite luminosity bins, and the Eddington ratio distribution is biased 
low by the same amount. A radio quasar subsample (with 1.5 < z < 2.3) has mean virial BH mass larger by 
~ 0.12 dex than the radio-quiet sample matched in luminosity and redshift. A broad absorption line (BAL) 
quasar subsample (with 1.7 < z < 2.2) has a virial mass distribution identical to that of the non-BAL quasar 
sample matched in luminosity and redshift, with no mean offset. 

Subject headings: black hole physics - galaxies: active - galaxies: fundamental parameters - galaxies: high- 
redshift - quasars: general - surveys 



1. INTRODUCTION 

There has been increasing interest in recent years in the role 
that supermassive black holes (SMBHs) play in galaxy forma- 
tion, primarily because of the discovery that most, if not all, 
present-day massive galaxies harbor a SMBH in their nuclei 
(e.g., Kormendy & Richstone 1995; Richstone et al. 1998), 
and that the mass of the nuclear SMBH is related to the bulge 
mass/luminosity (e.g., Magorrian et al. 1998) and even more 
tightly to the bulge stellar velocity dispersion (e.g., Ferrarese 
& Merritt 2000; Gebhardt et al. 2000a; Tremaine et al. 2002). 
These tight correlations imply that the formation of galaxies 
and the growth of the central SMBHs are ultimately connected 
(e.g.. Silk & Rees 1998; Kauffmann & Haehnelt 2000; Wyithe 
& Loeb 2003; Di Matteo et al. 2005; Hopkins et al. 2006). 

On the other hand, it has long been suggested that active 
galactic nuclei (AGN) or quasars are SMBHs in the process 
of accretion and growth (e.g., Salpeter 1964; Zel'dovich & 
Novikov 1964; Lynden-Bell 1969). Indeed, the integrated 
luminosity density of optically-selected AGNs, which repre- 
sents the accretion history of black holes, is consistent with 
the mass density in the local dormant SMBH population (the 
Soltan [1982] argument; e.g., Salucci et al. 1999; Yu & 
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Tremaine 2002). This is particularly encouraging, since with 
ever larger AGN/quasar samples in modern surveys, we can 
hope to understand the cosmic evolution of SMBHs within 
the framework of hierarchical structure formation, and shed 
light on the coevolution of SMBHs and their host galaxies. A 
central issue in this regard is to measure the masses of both 
inactive SMBHs and active AGNs/quasars. 

Perhaps the most reliable way to measure the black hole 
mass is via stellar/gas dynamics in the vicinity of the SMBH 
where its gravity dominates the dynamics (e.g., Richstone et 
al. 1998), but the relevant scales can be resolved only for the 
nearest galaxies. A more indirect mass indicator uses the ob- 
served tight correlations between the SMBH mass Mbh and 
the stellar bulge velocity dispersion a (the M-a relation), 
or the bulge luminosity (the M-L relation). Accurate mea- 
surements of the bulge velocity dispersion or bulge luminos- 
ity are difficult in luminous quasars, where the AGN light 
overwhelms that of the host galaxy; these relations are mostly 
applied to measure the local dormant SMBH mass function. 
More seriously, the M-a and M-L relations have not been 
measured directly at the most massive end, and the high-mass 
end of the BH mass functions estimated using the M-a and 
M-L relations differ by as much as one order of magnitude 
(Tundo et al. 2007; Lauer et al. 2007). Thus there are still 
unsettled issues on the usage of these techniques. 

A third method, reliant on AGN physics rather than galaxy 
properties, is reverberation mapping, which uses the tempo- 
ral offset between continuum and emission line variability to 
determine the distance R from the central engine to the broad 
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emission line region (BLR; e.g., Blandford & McKee 1982; 
Peterson 1993; Kaspi et al. 2000; Peterson et al. 2004). Us- 
ing the observed line width V and assuming that the BLR 
is virialized allows determination of the black hole mass: 
Mbh — G~^RV^. In a handful of cases, measurements of V and 
R for several lines in the same object have shown V oc 
consistent with the virial hypothesis (e.g., Peterson & Wandel 
2000). Furthermore, for the few cases in which we have more 
than one mass indicator, the reverberation mapping mass is 
also consistent with the dynamical mass or that derived from 
the Mbh - cr relation (e.g., Gebhardt et al. 2000b; Ferrarese et 
al. 2001; Nelson et al. 2004; Onken et al. 2004, 2007; Greene 
& Ho 2006; Davies et al. 2006). The reverberation mapping 
method is very time consuming, and we have reliable masses 
with this method for only about three dozen AGNs. 

However, reverberation mapping data has revealed a cor- 
relation between BLR size R and luminosity (Kaspi et al. 
2000, 2005), which allows us to estimate BH masses based on 
single-epoch spectra. Using this R-L correlation and rever- 
beration mapping masses, various empirical scaling relations 
(which we term virial BH mass estimators from now on) have 
been derived using H/3 (Kaspi et al. 2000; Vestergaard 2002; 
McLure & Jarvis 2002; Vestergaard & Peterson 2006), Ha 
(Greene & Ho 2005), Mgii (McLure & Jarvis 2002; McLure 
& Dunlop 2004), and Civ (e.g., Vestergaard 2002; Vester- 
gaard & Peterson 2006). This is the only practical method for 
measuring black hole masses for large numbers of objects at 
high redshifts and luminosities (where variability timescales 
are too long to measure easily; see the discussion in Kaspi 
et al. 2007), but the R-L relation has only been established 
in the luminosity range ^ 

1042 _ 1046 erg s-i and for redshifts 

z<0.3. 

The virial BH mass estimators are usually expressed as: 

log =a+blos — J-. +21og 7- 

^\ Mq J ^V1044ergs-V ^Vkms-i 

(1) 

where XLx, the surrogate for the BLR size R, is the continuum 
luminosity near the emission line (5100A for H/3, 3000A for 
Mgii and 1350A for Civ), and FWHM is the full-width-at- 
half-maximum of the line. Some authors prefer to use other 
quantities in their virial estimates. For example, Wu et al. 
(2004) suggest using recombination line luminosities instead 
of continuum luminosities because the latter might be con- 
taminated by jet emission or host galaxy starlight (also see 
Greene & Ho 2005). Peterson et al. (2004) and Collin et al. 
(2006) argue that the second moment of the line profile criine is 
a better quantity to characterize the emission-line widths than 
the FWHM. As we will show below, different measurement 
techniques yield systematic differences in derived quantities, 
and so it is important to use the original definitions of line 
widths and luminosities for whichever calibration is used. 

The UV Civ virial mass estimator can be used for z > 2 
quasars, as the Civ line is the only relevant line available 
in optical spectra in this redshift range. However, it has 
been suggested by a few authors that the Civ line might be 
a less secure mass estimator than Mgii or H/?. The Civ line 
tends to be asymmetric and blueshifted with respect to lower- 
ionization lines such as H/3 or Mgii (e.g., Gaskell 1982; Tytler 
& Fan 1992), and the most blueshifted objects tend to have the 
largest FWHMs (Richards et al. 2002b). These features sug- 
gest that Civ might be more severely affected than the other 
lines by a non-virialized gas component, and the measured 
line width could depend on the viewing angle, biasing the 



BH mass estimates. Indeed, using a sample of 80 z < 0.5 
Palomar-Green (Green, Schmidt, & Liebert 1986) quasars, 
Baskin & Laor (2005) showed that the Civ FWHM is poorly 
correlated with the H/3 FWHM (e.g., see their figure 3), sug- 
gesting that the two lines have different origins. On the other 
hand, the Mgii FWHM is well-correlated with that of H/3, and 
these two low-ionization line estimators usually give consis- 
tent virial masses (e.g., McLure & Dunlop 2004; Salviander et 
al. 2007). One of the purposes of this paper is to explore how 
the Civ estimator might be biased relative to low-ionization 
line estimators such as Mgii using a large sample of quasars 
from the Sloan Digital Sky Survey (SDSS; York et al. 2000). 

Despite the caveats in these virial estimators, virial BH 
masses have been measured for various AGN/quasar sam- 
ples covering a wide range of redshifts and luminosities. 
For example, McLure & Dunlop (2004) measured virial BH 
masses for 12698 quasars with 0.1 < z < 2.1 from the SDSS 
DRl quasar catalog (Schneider et al. 2003) using H/3 and 
Mgii. Vestergaard (2004) used H/3 and Civ to measure virial 
BH masses for a hybrid sample including 87 z < 0.5 Bright 
Quasar Survey (BQS, Schmidt & Green 1983) quasars, 114 
1-5 ^ z < 3.5 quasars, and ^ 150 z > 3.5 SDSS quasars. 
Greene & Ho (2007) present the BH mass function for ^ 8500 
broad-line AGNs with z < 0.35 from the SDSS. Fine et al. 
(2006) used composite spectra to measure the redshift evolu- 
tion of the mean BH mass for the 2QZ quasar sample (Croom 
et al. 2004) from z 0.5 to z 2.5, and Kollmeier et al. 
(2006) presented virial BH masses for a sample of 407 AGN 
with z ^ 0.3-4 selected from the AGN and Galaxy Evolu- 
tion Survey (AGES; Kochanek et al. 2004) and focused on 
the distribution of Eddington ratios. Yet by far the largest 
quasar sample is the recently published SDSS DR5 quasar 
catalog (Schneider et al. 2007), containing more than 77,000 
quasars, about half of which are homogeneously selected 
(e.g., Richards et al. 2006a; Shen et al. 2007). In this paper we 
present virial BH mass estimates for this optical quasar sam- 
ple and explore possible biases in their estimation. Though 
the dynamical range in luminosity and color of SDSS quasars 
is limited, the large size of the sample provides unprecedent 
statistics. 

The paper is organized as follows: in §|2]we describe our 
quasar sample. We present our procedures of estimating virial 
BH masses in § [3] where we compare results from different 
lines in detail. The distributions of black hole masses and Ed- 
dington ratios are described in § H] and we discuss our results 
as well as some general issues with virial estimators in ^ We 
summarize our results in §|6] Throughout this paper we adopt 
a flat ACDM cosmology: flM = 0.26, = 0.74 and h = Q.l\ 
(Spergel et al. 2007). 

2. THE SAMPLE 

The SDSS uses a dedicated 2.5-m wide-field telescope 
(Gunn et al. 2006) with a drift-scan camera with 30 2048 x 
2048 CCDs (Gunn et al. 1998) to image the sky in five broad 
bands {ugriz; Fukugita et al. 1996). The imaging data are 
taken on dark photometric nights of good seeing (Hogg et 
al. 2001), are calibrated photometrically (Smith et al. 2002; 
Ivezic et al. 2004; Tucker et al. 2006) and astrometrically (Pier 
et al. 2003), and object parameters are measured (Lupton et al. 
2001; Stoughton et al. 2002). Quasar candidates (Richards et 
al. 2002a) for follow-up spectroscopy are selected from the 
imaging data using their colors, and are arranged in spec- 
troscopic plates (Blanton et al. 2003) to be observed with a 
pair of double spectrographs. The quasars observed through 
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the Third Data Release (Abazajian et al. 2005) have been 
cataloged by Schneider et al. (2005), while Schneider et al. 
(2007) extend this catalog to the DR5 (Adelman-McCaithy et 
al. 2007a). 

Our parent sample is the published SDSS DR5 quasar cat- 
alog (Schneider et al. 2007), which contains 77,429 bona fide 
quasars that have luminosities larger than M,- = -22.0 (us- 
ing a slightly different cosmology in that paper) and have 
at least one broad emission line (FWHM > 1000 km s"') or 
have interesting/complex absorption features. About half of 
the quasars in this catalog were selected from a uniform al- 
gorithm (as described in Richards et al. 2002a), which is flux 
limited to^ / = 19.1 at z < 3 and / = 20.2 at z > 3 (magni- 
tudes are corrected for Galactic extinction using the Schlegel, 
Finkbeiner & Davis 1998 map). These objects are flagged 
with the UNIFORM FLAG = 1 in the DR5 quasar catalog, 
and they can be used to construct a statistically homogeneous 
sample (e.g., Richards et al. 2006a; Shen et al. 2007). 

There are also flags indicating whether or not a quasar is 
detected in the FIRST (Becker et al. 1995) or ROSAT (Voges 
et al. 1999) surveys, which we will use to define our radio 
quasar subsample (see 34.6b . Most of our analysis will avoid 
broad absorption line quasars (BALs), whose line widths are 
problematic to measure; we will come back to the BALs in 
SM We have identified - 4200 BALs (with 1.7 < z < 4.2 
for Civ and 0.5 < z < 1.9 for Mgii; Shen et al. 2008) in the 
DR5 quasar catalog using traditional BAL criteria (Weymann 
et al. 1991). This list of BALs is by no means complete, and 
an official DR5 BAL catalog is forthcoming (Hall et al., in 
preparation; see Trump et al. 2006 for the DR3 BAL catalog). 

All spectra are reduced with the new version of the SDSS 
spectroscopic reduction pipeline, as described in the DR6 pa- 
per (Adelman-McCarthy et al. 2007b); the flux scale of these 
spectra is higher than that of previous releases by roughly 
38%. Thus in this paper all luminosities measured from spec- 
tra will be typically larger by ^ 0. 14 dex than previous values 
for the same objects. 

A comprehensive catalog of our measurements of the spec- 
tral properties of each quasar is given in the electronic version 
of this paper; the columns are described in Table 1. This cat- 
alog can be regarded as an extension to the published DR5 
quasar catalog. 

3. BLACK HOLE MASS ESTIMATES 

To estimate BH masses using virial estimators one needs 
two parameters: the width of an emission line and the corre- 
sponding continuum luminosity. However, different authors 
have used different definitions of line width in their calibra- 
tions of the virial mass estimators. In what follows, we take 
care to define line width in a way consistent with each cali- 
bration we use. We also note that the assumed virial coeffi- 
cient /, which accounts for our ignorance of the BLR geom- 
etry, is different in different calibrations. For instance, / = 1 
in McLure & Dunlop (2004) and / ~ 5.5 (e.g., Onken et al. 
2004) in Vestergaard & Peterson (2006). Thus different ver- 
sions of calibration for the same line do not necessarily yield 
the same results, especially when extrapolated to high lumi- 
nosity quasars. We will discuss this further in 35.11 

We focus on nonBALs throughout this section and 34. U 

^ There are a few i > 19.1 quasars at z < 3 which were selected by 
the high-z (griz) branch of the targeting algorithm (Richards et al. 2002b). 
The fraction of these objects is tiny (< 2%). Although these objects have 
UNIFORM FLAG = 1 in Tab le 1, they will be rejected when constructing a 
flux-limited sample (see ^4.3t . 
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Fig. 1 . — Examples of spectral fits. From top to bottom: H/3, Mgll, CIV. 
In each panel, the upper and lower black solid lines shows the original and 
continuum+Fe subtracted spectra; the gray lines show the flux density errors; 
the red lines show the full fits; the orange lines show the fitted power-law 
continuum; the lower blue lines in the upper and middle panels show the iron 
template fits; the magenta lines show the fits for the emission lines and green 
lines show the Gaussian components. The blue hoiizontal segment in the 
bottom panel marks the FWHM for CIV. 

I4.5l because of the ambiguities of line wid th m easurement for 
BALs; we return to the issue of BALs in 34.61 

3.1. Line widths and continuum luminosities 
3.1.1. HPandMgU 

For the H/3 and Mgii estimators, we have adopted the 
calibrations in McLure & Dunlop (2004), hence a = 0.672, 

= 0.61 for H/3 and a = 0.505, b = 0.62 for Mgii in equation 
([T]). The H/3 estimator has been calibrated directly with rever- 
beration mapping, while the Rblr - relation for Mgii at 
3000A is empirically determined in McLure & Jarvis (2002) 
based on H/3 reverberation mapping masses. We note that the 
calibration for Mgii in McLure & Dunlop (2004) is an up- 
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TABLE 1 
Catalog Format 



Column Format Desciiption 



1 


A18 


^T^S^ ril? S Hp*;! onntinn hhmTTi';'; cc-l-ririmmcc f^y^OOO 0^ 


2 


F11.6 


Right ascension in decimal degrees (J2000.0) 


3 


F11.6 


Declination in decimal degrees (12000. 0) 


4 


F7.4 


Redshift 


5 


F8.3 


Mi(z = 2}(h = 0.71, Um = 0.26, = 0.74, /T-coiTected to j = 2, following Richards et al. 2006a) 


6 


F7.3 


Bolometric luminosity (log(L(,oi/ei'gs"')) 


7 


15 


Spectroscopic plate number 


g 


15 


Spectroscopic fiber number 


9 


16 


MJD of spectroscopic observation 


10 


112 


Target selection flag when the spectrum was taken (i.e., using TARGET photometry) 


11 


13 


F/R5r selection flag (0 or 1) 


12 


13 


ROSAT selection flag (0 or 1) 


13 


13 


Uniform selection flag (0 or 1) 


14 

14 


Ij 


i5/\ij nag [y) or i j 


15 


17 


H/3FWHM(kms-') 


16 


F8.3 


Monochromatic luminosity \Lx at 5100A (lO'*'' erg s"') 


17 


F7.3 


Virial BH mass estimated using H/3 (log(MBH.vir/'W0)) 


18 


17 


MgIIFWHM(kms-') 


19 


F8.3 


Monochromatic luminosity \L\ at 3000 A (lO'*'* erg s"') 


20 


F7.3 


Virial BH mass estimated using Mgll (log(MBH,vii7M© )) 


21 


17 


ClVFWHM(kms-') 


22 


F8.3 


Monochi'omatic luminosity \Lx at 1350A (lO**** erg s"') 


23 


F7.3 


Virial BH mass estimated using CIV (log(MBH,vir/W0)) 


24 


F7.3 


Virial BH mass (using H/3 for z < 0.7; Mgll for 0.7 < z < 1.9 and CIV for z > 1.9) 


25 


17 


CIV-Mgll blueshift (km s"') 


26 


F8.3 


Mean spectrum S/N (signal-to-noise ratio) 



Note. — (1) Objects in this table are in the same order as in the DR5 quasar catalog (Schneider et al. 2007); (2) A"— corrections 
are the same as in Richards et al. (2006a); (3) Bolometric luminosities are computed using bolometric corrections in Richards 
et al. (2006b) using one of the 5100A, 3000A, or 1350A monochromatic luminosities depending on redshift; (4) Entries reading 
-9999 for FWHM, luminosity, BH mass or CIV-Mgll blueshift measurement indicate that this quantity was not measureable 
from the SDSS spectrum, either because it fell outside of the SDSS spectral coverage or because of low S/N. 



dated version of that in McLure & Jarvis (2002), where in the 
former only objects with 10^"^ erg s~' < \L\ < 10'^^ erg s~' 
were included in their fitting. This luminosity range is suit- 
able for the SDSS quasar sample studied here. 

We follow McLure & Dunlop (2004, Appendix A) to mea- 
sure the line FWHM for H/3 and Mgii with slight changes in 
details: 

1) a power-law continuum and an iron emission line tem- 
plate are simultaneously fitted to the spectrum with the 
H/3 or Mgii emission line regions excluded. The fitting 
ranges are [4435,4700] A and [5 100,5535] A for H/3, 
[2200,2675]A and [2925,3090]A for Mgii. The iron 
template for H/3 is taken from Boroson & Green (1992) 
and the iron template for Mgii is a modified version of 
the Vestergaard & Wilkes (2001) template (Salviander 
et al. 2007) which extends under the Mgii line itself. In 
the fitting, the normalization and velocity broadening 
of the iron template are left as free parameters. 

2) The best-fit continuum and iron emission are then sub- 
tracted from the spectrum, and two Gaussians are fit- 
ted to each emission line, one for the broad line com- 
ponent and the other for the narrow line component. 
The line fitting ranges are [4700,5 100] A for H/3 and 
[2700,2900]A for Mgii. The FWHM of the narrow 
component is constrained to be less than 1200 km s"', 
and that of the broad component is constrained to be 
larger than 1200 km s"' (e.g., Hao et al. 2005). In the 
case of H/3, two additional Gaussians whose FWHMs 
are tied to that of the narrow H/3 component are fit- 
ted simultaneously for [Oiii] 4959A and [Oiii] 5007A. 



The FWHM of the broad Gaussian-component of H/3 
and Mgii is then taken as the line width to be inserted 
in equation ([T]i. 

It is not entirely clear if there exists a strong narrow line 
component for Mgii, thus the way to measure the Mgii 
FWHM in McLure & Dunlop (2004) and here is not justi- 
fied. If we measure the Mgii FWHM from the two-Gaussian 
fitted flux, it will typically be smaller by 0.15 dex. How- 
ever, in order to use their calibration we have to follow the 
same procedure to measure the Mgii FWHM. As we will see 
later in ^3.3.11 the fitted broad line FWHMs for Mgii and H/3 
are quite similar, and as a result, both line estimators yield 
consistent virial masses. 

For continuum luminosities \L\ at 5100A for H/3 and at 
3000A for Mgii, we have used the fitted continuum flux den- 
sity at the corresponding wavelength and corrected for Galac- 
tic extinction using the Schlegel, Finkbeiner & Davis (1998) 
map. 

The calibration of the H/3 mass estimator is not unique, 
and other authors have determined slightly different versions 
(e.g., Vestergaard 2002; Vestergaard & Peterson 2006). We 
find other forms of the H/3 estimator yield comparable but 
not identical results. In particular, the Vestergaard & Peterson 
(2006) H/3 calibration gives a systematic ~ 0.15 dex offset 
compared with the McLure & Dunlop (2004) calibration. We 
will come back to this point in |5] 

3.L2. Civ 

For the Civ estimator, we use the calibration in Vestergaard 
& Peterson (2006), which has a = 0.66 and b = 0.53 in equa- 
tion ([B. 
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2 3 

Redshift 

Fig. 2. — Left: correspondence between M,(z = 2) and Lboi- Lines show median values and standard deviations. Tlie absolute magnitude has been K-corrected 
to the i'-band measured at z = 2 (i.e., a rest-frame wavelength of 2500 A), following Richards et al. (2006a). Right: redshift distribution of for our sample, 
where black points are uniformly selected quasars while gray points are quasars selected with special targeting algorithms. 




Fig. 3. — Distiibutions of FWHMs as functions of redshift and luminosity. The FWHMs for the three lines are plotted together. Contours are local point 
density contours, drawn to guide the eye. In the left panel, we also show the mean and 1 —cr of the FWHM distribution as a function of redshift. FWHMs show 
little or no dependence on either redshift or luminosity, and their distributions follow a log-normal at all redshifts and luminosities (see Fig.|4). 



Iron emission contamination is not a problem for Civ, so 
we have simply fitted a power law continuum to regions in 
the vicinity of Civ, namely [1445,1465]A and [1700,1705]A 
. This continuum is then subtracted from the spectrum for the 
line fitting which follows. 

The FWHM for Civ as defined in Vestergaard & Peterson 
(2006) is indeed the "full-width-at-half maximum" of the full 
line profile, whether or not the line profile is Gaussian. Since 
many Civ lines have profiles "peakier" than a Gaussian, the 
FWHM obtained from a single Gaussian fit overestimates its 
formal definition. For this reason, following the standard ap- 
proach in the literature (e.g., Baskin & Laor 2005; Fine et al. 
2006), we use a set of three Gaussians to fit the Civ line re- 



gion [1500,1600]A. The FWHM is then measured from the 
fitted model flux, providing better measurements of FWHM 
for noisy spectra. The Civ line in high redshift quasars may 
be biased by associated absorption lines. To reduce their ef- 
fect, a second line fit is performed which excludes pixels that 
are 1 .5a below the first fitted model, where a is the estimated 
error in flux density. The model is replaced by the second fit 
if the reduced is smaller than in the first fit. Also, if one or 
two of the three Gaussians have integrated flux less than 5% 
of the total model flux, they are removed from the model flux; 
this step ensures that there are no unphysically narrow peaks 
to bias the FWHM measurement. As in the cases of H/3 and 
Mgii, we require FWHM > 1200 km s"' . 
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1-2 3-4 3-6 3-8 4.0 4.2 3-2 3.4 3.6 3.8 4.0 4.2 3.2 3.4 3.5 3.8 4.0 4.2 

0.1<z<0.4 0.4<z<0.7 0.7<z<1.1 

Fig. 4. — Distributions of log FWHM (in km s"' ) in different redshift (liorizontally arranged) and luminosity (vertically arranged) bins for H/3, Mgll and CIV. 
Solid histograms show the results without a S/N cut and dashed histograms show the results with S/N > 10. We only plot a histogram if there are more than 
200 objects (no S/N cut) or 50 objects (with S/N cut) in each bin. The number of quasars drops rapidly when we impose the S/N cut, but the solid and dashed 
histograms show almost identical distributions. The means and dispersions of a fitted log-normal to the distributions are shown in the upper-left comers of each 
panel. 



The Civ FWHMs measured in this way are systematically 
narrower by ^ 0.15 dex than what results from a single Gaus- 
sian fit. Clearly, we would get systematic offsets in the Civ 
based BH masses if we used the single Gaussian FWHM. 

The 1350A continuum luminosity XL\ is measured us- 
ing the mean flux density within the wavelength range 
[1340, 1360] A, and corrected for Galactic extinction. Given 
the wavelength coverage of the SDSS spectra, this implies 
that we cannot estimate Civ black hole masses for quasars 
with redshift below 1.8. 

3.2. Bolometric Luminosities 
For each object in our sample that has a measurable virial 
mass, we use its continuum luminosity XL\ (at 1350A for 
Civ, 3000A for Mgii and 5 lOOA for H/3) to estimate the bolo- 
metric luminosity. The bolometric correction factors (BC) are 
computed using the composite SED for a sample of SDSS- 
DR3 quasars constructed by Richards et al. (2006b), and are 
slightly different from the commonly adopted values in Elvis 
et al. (1994). In particular, BC1350 = 3.81, BC3000 = 5.15, 
BC5100 = 9.26. When two continuum luminosities are avail- 
able, they give quite similar bolometric luminosities. The typ- 
ical error in bolometric corrections assuming a universal SED 
is 0.1 dex for optically-selected SDSS quasars (Richards et 
al. 2006b), but it might be sHghdy larger for BCbsq. The left 
panel of Fig. |2] shows the relation between Lboi and the z = 2 
/-band absolute magnitude calculated using the /T-correction 
in Richards et al. (2006a), while the right panel shows Lboi 
versus redshift for our sample (the darker points are the uni- 
formly selected subsample). 



Fig. [3] shows the distribution of log FWHM in bins of red- 
shift and bolometric luminosity. The FWHMs of all three 
emission lines are only weakly dependent on either redshift or 
luminosity within our sample, and they follow a log-normal 
distribution with typical dispersion ~ 0.1-0.2 dex. This 
width is real and is not dominated by noise; we find essen- 
tially identical distributions when restricting our sample to 
objects with S/N per pixel greater than 10, or greater than 20. 
To further demonstrate that the FWHM distributions do not 
evolve. Fig. |4] shows the FWHM distributions for the three 
lines in different redshift and luminosity bins. Our measured 
distributions of FWHMs are consistent with other investiga- 
tions (e.g., Baskin & Laor 2005; Salviander et al. 2007), with 
similar widths and a '-^ 0.1 dex difference in the mean values^ 
(which might be caused by small differences in the techniques 
to determine FWHMs). The observed distribution of FWHMs 
has profound effects on our determination of the distribution 
of virial BH masses, as will be further explored in section |43] 

3.3. Comparison of Emission Lines 

There are subsamples of quasars for which we have two 
lines available; given the SDSS spectroscopic coverage, we 
observe both H/? and Mgii for 0.4 < z < 0.9, and both Mgii 
and Civ for 1 .7 < z < 2.2. Therefore we can check if different 
BH mass estimators give consistent results. 

3.3.1. HP vs Mgll 



The difference of 0. 1 dex in FWHM will result in ~ 0.2 dex mean 
offset in virial masses. We do not consider this substantial, since the zero- 
point of the virial calibrations has uncertainties of the same level. 
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Fig. 5. — Comparisons between the H/3 and the Mgll estimators, left: histogram of the FWHM ratios, right: histogram of the virial BH mass ratios. 



10.0 



9.5 



9.0 



rS 8.5 



cn 
o 




o 




7.5 - 



7.0 



8.5 



7.0 



7.5 



8.5 



9.0 



9.5 10.0 



8.5 



Log M^S (Mo) 



9.0 9.5 

Log (Mo) 



10.0 



Fig. 6. — A one-to-one comparison between two virial mass estimators, left: H/3 versus Mgll; right: Mgll versus CIV. The congelation is much better between 
H/3 and Mgll than the one between Mgll and CIV. Contours are local point density contours, drawn to guide the eye. 



For a subsample of ^ 8,000 quasars with redshift 0.4 < 
z < 0.9 we have both H/3 and Mgii Hnewidth measurements. 
The distribution of log(H/3 FWHM/Mgii FWHM) is shown 
in the left panel of Fig. |5] which follows a Gaussian with 
mean 0.0062 and dispersion 0.11 dex. The distribution of 
log(Mg^/M^]|'') is shown in the right panel of Fig. |5] which 
follows a Gaussian with mean 0.034 and dispersion 0.22 dex. 
The small discrepancy in the mean offsets between our re- 
sults and those in McLure & Dunlop (2004, where the mean 
FWHM ratio and BH mass ratio offsets ai-e -0.004 and 0.013 
dex, respectively) is perhaps due to different samples, but is 
totally negligible considering the ^ 0.3-0.4 dex uncertainties 
of either of the BH mass estimators. Therefore, the Mgii and 
H/3 estimators give consistent BH masses. 

A plot of the H/3 FWHM against the Mgii FWHM is very 
similar to Fig. 3 in Salviander et al. (2007); the two FWHMs 
are correlated, with a slope differing slightly from unity. 



We plot the H/J against the Mgii virial masses in the left- 
hand panel of Figure |6l The two mass estimators correlate 
strongly with one other over two orders of magnitude in mass. 

3.3.2. M^ii vs Civ 

Linewidths of both Mgii and Civ are available for a sub- 
sample of ^ 15,000 quasars with redshift 1-7 < z < 2.2, of 
which ^ 5, 000 quasars have both Mgii and Civ virial masses 
available. However, before we compare the results, we review 
some of the unique characteristics of the Civ line. 

First, the definitions of FWHM are quite different for 
Mgii and Civ. While the broad component of Mgii is well- 
described as a Gaussian in our fits, the vast majority of Civ 
lines are by no means Gaussian. The difference in line shapes 
suggests that the broad Mgii and Civ might have different 
physical origins. Second, it is well known that high ionization 
lines such as Civ often show systematic blueshifts with re- 
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Fig. 7. — Various properties against the CIV-Mgll blueshift. Contours are local point density contours, drawn to guide the eye. Left: CIV FWHM against 
blueshift. The mean value of FWHMs rises with increasing blueshift. For those objects with the largest blueshift, their FWHMs rarely drop below 4000 km s"' . 
Middle: Mgll FWHM against blueshift. In contrary to the CIV case, the Mgll FWHM shows a mild decrease with increasing blueshift, but it is also consistent 
with no correlation at all considering the huge scatter. Right: Bolometric luminosity against blueshift. The mean luminosity increases with blueshift, but again 
the dispersion is large. 

spect to low ionization lines such as Mgii (e.g., Gaskell 1982; 
Tytler & Fan 1992; Richards et al. 2002b). Richards et al. 
(2002b) used 800 SDSS quasars to show that Civ is sys- 
tematically blueshifted with respect to Mgii by ^ 900 km s"' 
on average. They interpreted this apparent blueshift as due 
to the absorption of the Civ red wing, as might be expected 
if the broad-line region gas is actually a radiatively driven, 
equatorial outflow. Such disk-wind models are compelling 
for a variety of reasons, particularly their success in explain- 
ing the properties of broad absorption line quasars (e.g., Mur- 
ray et al. 1995; Proga et al. 2000; Elvis 2000). Richai'ds et 
al. also see a trend of increasing average Civ FWHM with 
increasing Civ-Mgii blueshift, which again suggests that the 
broad Civ line is dominated by a disk wind component; in this 
case the Civ-Mgii blueshift reflects the orientation of the line- 
of-sight (LOS). There is additional evidence supporting this 
disk wind/orientation scheme, including the anti-correlation 
between the Civ-Mgii blueshift and radio loudness, and the 
correlation between the Civ-Mgii blueshift and BAL fraction 
(see also Gallagher et al. 2005). These properties are also seen 
in a sample of low redshift broad line AGNs (Sulentic et al. 
2007). Of course, alternative interpretations of these features 
are possible, which depend more on the structure of the wind 
than on the LOS orientation (see Richards 2006). 

If the broad Civ is indeed contaminated or even dominated 
by a disk wind component, then the measured Civ FWHM 
does not represent the actual virial velocity of the BLR. If the 
LOS is along the disk wind, then the FWHM will overesti- 
mate the virial velocity, while if the LOS is along the pole 
direction, the FWHM will probably underestimate the virial 
velocity. Although the Civ estimator is calibrated in such 
a way that it gives, on average, BH masses consistent with 
those estimated from H/3 (e.g., Vestergaard 2002; Vestergaard 
& Peterson 2006), there will be biases in the Civ-based virial 
masses if the above picture is correct. 

We explore these issues by plotting various quantities 
against the Civ-Mgii blueshift in Fig. |7] The left panel shows 
the coiTelation between the Civ FWHM and the Civ-Mgii 
blueshift. We see the rise of the average Civ FW^HM as the 
blueshift increases, as reported by Richards et al. (2002b). In 
particular, very few objects with large blueshifts have FWHM 
less than 4000 kms"', while the distribution of FWHM is 
broader for objects with small blueshifts. The middle panel 



shows a mild anti-correlation between the Mgii FWHM and 
the Civ-Mgii blueshift, and the right panel shows a mild 
correlation between the bolometric luminosity Lboi and the 
blueshift, which has been found before (e.g., Richards et 
al. 2002b). The cause of the anti-correlation between Mgii 
FWHM and Civ-Mgii blueshift is not clear at this point. It 
might reflect inclination effects as well. The different corre- 
lations of the Civ and Mgii line widths with blueshift suggest 
that the two lines have different origins, as we argued above. 

Since the Civ estimator is calibrated using H/3 reverbera- 
tion mapping masses (e.g., Vestergaard 2002; Vestergaard & 
Peterson 2006), there should be no mean offset in either the 
FWHMs or the estimated BH masses, as we find (see the black 
solid histograms in the bottom panels in Figs. [8]and|9]l. While 
the log of the ratio of Mgii to Civ FWHMs follows a Gaus- 
sian with mean 0.027 and dispersion 0.18 dex, the correlation 
between the two line widths is very weak (top panel of Fig. 
[SJ; indeed the width of 0.18 dex is larger than the dispersion 
in the Mgii FWHM distribution itself. Correspondingly, the 
log of the ratio of Mgii to Civ BH masses follows a Gaussian 
with mean -0.06 dex and dispersion 0.34 dex. Baskin & Laor 
(2005) similarly found weak correlations between the Civ and 
H/3 FWHMs at low redshift, while Netzer et al. (2007) report 
a similar effect at z w 2 (but see Kelly & Bechtold 2007 for 
an alternate view). Moreover, this scatter is systematic with 
blueshift. The blue points and histogram in Figure prefer to 
objects with blueshift larger than 1000 km s"', while the red 
points have smaller blueshift. 

The Civ and Mgii virial mass estimators are plotted against 
one another in the right-hand panel of Figure |6] The corre- 
lation is much weaker than we found for the Mgii-H/3 com- 
parison. The ratio of these two masses is correlated with the 
Civ-Mgii blueshift, as is shown in Fig. |9] (top panel). Given 
the similarity between the properties of the Mgii and H/3 lines, 
and the evidence that the Civ line may be affected by a disk 
wind, we argue that H/J and Mgii are better BH mass esti- 
mators than is Civ. In particular, we emphasize the danger 
of using Civ as a BH mass estimator for individual objects 
and small samples. For example, Kurk et al. (2007) and Jiang 
et al. (2007a) have obtained observations of both Mgii and 
Civ in four z ^ 6 quasars. In the three quasars in which Civ- 
Mgii blueshifts are observed, the Civ virial mass is ^ 2 times 
larger than the Mgii mass, while they are in good agreement 
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Fig. 8.— Top: CIV FWHMs against MgllFWHMs for a sample of quasars 
with botli FWHMs available. Blue and red dots show those with CiV-Mgll 
blueshift greater and less than 1000 km s"' ; the solid line shows the unity 
relation. Bottom: histograms of the ratios of the Mgll FWHM to the CIV 
FWHM. Black histogram is for the whole sample; blue and red ones are for 
objects with CIV-Mgll blueshift greater and less than 1000 km s"' . 

for the one object with a Civ-Mgii redshift, consistent with 
the trends we see in Figure |9] Finally we point out that the 
ratio of Mgii and Civ virial masses does not depend on lumi- 
nosity in the range where both line estimators are available, 
excluding luminosity-based biases. 

Of course we cannot measure the blueshift, and thereby cor- 
rect the Civ mass estimate for the bias of Figure|9] of objects 
for which our spectral coverage only includes Civ. However, 
Richards et al. (2002b) show that the Civ-Mgii blueshift cor- 
relates broadly with the Civ line asymmetry, suggesting that 
we might use this asymmetry as a surrogate to estimate the 
mass estimate bias. However, the correlation between the 
Civ-Mgii blueshift and the Civ line asymmetry shows a large 
amount of scatter, making it difficult to correct for this bias 
using the Civ line asymmetry. 
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Fig. 9. — Top: Difference in the Mgll and CIV based virial BH masses as 
function of the CIV-Mgll blueshift. For those with small (large) blueshifts, 
the CIV estimator systematically underestimates (overestimates) the BH 
mass as compared to the Mgll based virial masses. Contours are local point 
density contours, drawn to guide the eye. Bottom: Distributions of the ratios 
of the two virial mass estimates. The black histogram is for the entire sam- 
ple; blue and red ones are for objects with CIV-Mgll blueshift greater and 
less than 1000 km s"' . 

It is beyond the scope of this paper to carry out a thor- 
ough investigation of the BLR and disk geometry and dynam- 
ics needed to understand the physics behind the trends we've 
found. Therefore, in the following sections, we simply use the 
Civ estimator to measure BH masses at z > L9. Our confi- 
dence in the Civ estimator is built on the large sample studied 
here and the belief that there is no systematic offsets in the 
mean using the Civ estimator 

4. DISTRIBUTION OF BLACK HOLE MASSES 

We present virial BH masses for the SDSS DR5 quasars 
with z ^ 4.5. Excluding those objects for which we cannot ob- 
tain reliable FWHM or continuum luminosity measurements 
due to poor spectra or redshift limitations, we have 56,872 



10 



SHEN ET AL. 



objects (nonBALs) with measurable virial BH masses. The 
masses are based on the H/3 estimator for z < 0.7, Mgii for 
0.7 < z < 1.9, and Civ for 1.9 < z < 4.5, to avoid using mea- 
surements which are too close to the red and blue ends of the 
spectra. We should note that the quasars in this catalog were 
not all uniformly targeted; in particular, special targeting al- 
gorithms were used to select quasar candidates to fainter mag- 
nitudes than the main quasar survey (Richards et al. 2002a). 
The completeness of these fainter quasars is poor We start 
by neglecting this effect; we take all the quasars in the DR5 
quasar catalog with measurable masses and present the red- 
shift evolution of virial BH masses 04.1l i and Eddington ra- 
tio distribution s (34.2l l. We re strict ourselves to the complete 
subsample in ^4.31 and ^4.41 where we model the observed 
black hole mass distribution. We describe the active BH mass 
function in ^4.51 and discuss radio-loud and BAL quasars in 
SM 

4. 1 . Redshift evolution of virial BH masses 

Fig. [TOl shows the dependence of virial BH masses on red- 
shift, where gray dots represent individual objects and the 
filled circles show the mean value at each redshift. Be- 
cause of the incompleteness of the sample at the faint end, 
the "mean" BH masses here are weighted towards high lu- 
minosity objects at all redshifts. Nevertheless, the mean BH 
mass increases with redshift below z ~ 2, and then flattens 
out above z ^ 2; the most massive BHs seem to turn off 
around this epoch. There is clear evidence of a limiting BH 
mass ^ 3 X lO''- 10'" Mq, consistent with the conclusions 
of McLure & Dunlop (2004) for z < 2.1 SDSS DRl quasars. 
The trend of this figure is quite similar to fig. 4 of Vestergaard 
(2004), except that we have filled the gaps in their figure with 
our much larger sample. 

4.2. Distribution of Eddington ratios 

The relation between virial BH mass and bolometric lu- 
minosity is shown in Fig. [TT] The figure shows objects us- 
ing the three virial estimators in different colors (red for H/3, 
green for Mgii and blue for Civ), while the diagonal lines in- 
dicate various Eddington ratios L/L-Eid, where LEdd = 1.3 x 
10-'^(Mbh/A^o) erg s"' is the Eddington luminosity. The vast 
majority of quasars in our sample are bounded by Lgdd and 
O.OlLEdd, with a geometric mean value of ^ O.lLEdd- We are 
not surprised to see that objects using these virial relations 
occupy a stripe in the mass-luminosity plane, with a slope 

b ^ 0.5 and scatter ^ 0.2-0.3 dex around the mean. This 
simply reflects the fact that the FWHM distribution is almost 
independent on luminosity, and has a very narrow width of 
^ 0.1-0.15 dex at fixed luminosity (see Fig. |4|i. This also 
explains the remarkably similar appearance of the MfiH-i^boi 
diagram in several other investigations even though the un- 
derlying samples are quite different (e.g.. Woo & Urry 2002; 
Kollmeier et al. 2006), simply because similar virial estima- 
tors are used to estimate BH masses, and the FWHM distribu- 
tions of different samples are similar 

Fig. [T2I shows the distributions of Lboi/i^Edd in five lumi- 
nosity and four redshift bins. The distribution of Eddington 
ratios is log-normal in almost all of the bins. No appreciable 
redshift evolution is seen for objects with the same luminos- 
ity. The widths of the distributions are as small as ^ 0.3 dex, 
similar to those in Kollmeier et al. (2006). Since virial mass 
is approximately proportional to the square root of luminos- 
ity, the widths of the observed Eddington ratio and the virial 



BH mass distributions in each luminosity/redshift bin should 
be comparable, as is seen. The width of the virial BH mass 
distributions within each luminosity/redshift bin is thus ~ 0.3 
dex or less. 

At face value, the very narrow widths of the virial BH mass 
distribution in each luminosity bin seem somewhat surprising, 
since many studies have quoted significantly larger scatter 
between single-epoch virial masses and both reverberation- 
mapped masses on the one hand and stellar velocity disper- 
sions on the other. For instance, McLure & Jarvis (2002) 
find an rms scatter of 0.4 dex of virial Mgii masses around 
reverberation-mapping masses, while Vestergaard & Peterson 
(2006) find scatters of 0.3 and 0.5 dex of the Civ and H/3 
virial masses about the RM masses respectively. An indepen- 
dent estimate may be derived by comparing virial masses with 
stellar velocity dispersions; Greene & Ho (2006) find an rms 
scatter of 0.4 dex of virial Ha masses around the expectation 
from the M-a relation. In all of these cases, it is important 
to remember that the objects in question are typically much 
less luminous (and at much lower redshift) than the quasars 
we study here, but these are the measurements at our disposal 
at the present time. 

However, the observed width of the virial mass distribution 
within a given luminosity/redshift bin is not the uncertainty 
in the virial mass estimators themselves. Although we do not 
have a deterministic model for the origin of the uncertainty in 
the virial relations, it must result both from the imperfect rela- 
tion between luminosity and BLR radius on the one hand, and 
from the imperfect relation between FWHM and BLR virial 
velocity on the other. Thus the uncertainty in the virial mass 
estimators comes from the uncorrelated rms scatters in both 
luminosity and FWHM. When we impose luminosity cuts to 
our sample, i.e., either by working within luminosity bins or 
because of the flux limits of our sample, we will artificially 
narrow the observed distribution of virial masses. 

At the same time, if the underlying BH mass distribution is 
such that there are more BHs towards lower masses/lower lu- 
minosities, then more objects will scatter from low-mass bins 
to high than the other way around, giving rise to a Malmquist- 
type bias (e.g., Eddington 1913; Malmquist 1922) in the ob- 
served virial black hole mass and Eddington ratio. In the next 
section, we use Monte Carlo simulations to investigate the im- 
pact of luminosity cuts and Malmquist bias on our observed 
distributions of luminosity, line width, virial BH masses, and 
Eddington ratios. The results are shown in § 14.41 where we 
also present an analytic derivation of the Malmquist bias. 



4.3. Monte Carlo simulations 

We ask whether we can reproduce the observed distribu- 
tions given a set of reasonable assumptions for the true under- 
lying distributions of black hole masses and Eddington ratios. 
In particular, our model should be able to reproduce the ob- 
served distributions of FWHMs in each redshift-luminosity 
bin, and a reasonable amount of uncertainty in virial rela- 
tions. Also, the model should produce a luminosity distri- 
bution consistent with the observed quasar luminosity func- 
tion. Since the BH mass and Eddington ratio distributions 
may evolve with redshift, we consider two redshift ranges: 
0.7 < z < 1.0, where the Mgii estimator is used throughout, 
and 1.9 < z < 2.1, where the Civ estimator is used. We only 
consider the subset of uniformly-selected quasars in this sec- 
tion, and we restrict ourselves to objects with log ( ^'"l, ) > 
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Fig. 10. — The cosmic evolution of virial BH masses for our entire DR5 quasar catalog with measurable virial masses. The two dashed vertical lines mark the 
transitions from one estimator to another. The filled circles show the mean value in bins of redshift of width Az = 0.1 for z < 1.9 and Az = 0.2 at higher redshit. 
The gaps around z ~ 2.8 and z ~ 3.5 are caused by the inefficiency of selecting quasars around these redshifts (e.g., Richards et al. 2002a). The upper envelope 
flattens out at z > 2 but decreases towards lower redshift, which suggests that the most massive SMBHs are no longer shining as optical quasars at lower redshift. 
There is clearly an upper limit at all redshifts of black hole masses between 3 X lO' and lO'" Mq, shown as the two dash-dotted horizontal lines. 



45.9 and log (^^^pr^ > 46.6, respectively, above which our 

sample is close to complete in our two redshift ranges (e.g., 
Richards et al. 2006a). 

In outline, we begin by specifying a model for the true black 
hole mass distribution. For a given mass, we will assume a 
single underlying central Eddington ratio (and thus bolomet- 
ric luminosity) and a single underlying central FWHM. To 
get observed quantities (continuum luminosities and FWHM), 
we add scatter independently to the central Eddington ratio 
and FWHM. Our model is thus characterized by the assumed 
black hole mass function, the assumed relation between un- 
derlying Eddington ratios and black hole mass, and the scat- 
ters which take us to the observed luminosities and FWHM. 
We will present a model with reasonable (although not nec- 
essarily unique) choices for these parameters, which fits the 
data quite well. 

The observed quasar luminosity function has a power-law 
form at the bright end (Richards et al. 2006a; Hopkins et al. 
2007). We therefore start by assuming a power-law black hole 
mass distribution. For a given true black hole mass MBH.true, 
we assume a central underlying Eddington ratio: 



we can determine a central continuum luminosity. Given this 
luminosity and the true black hole mass, we assume that the 
virial estimator, equation with the observed values of a 
and b appropriate for the sample we're simulating, holds ex- 
actly for the central value of FWHM. That is, this central 
FWHM is given by: 



/FWHM\ 

A^BH.true 



log 



b{\og 



1044 



ergs 



(3) 



log 



= Ci+C2log 



BH.true 



(2) 



where Ci and C2 are constants to be set by the model. Given 
this Eddington ratio and our fiducial bolometric correction. 



This does not give us observed quantities yet. To go from 
the central to observed bolometric and continuum luminosi- 
ties, we add a scatter to each simulated object selected from 
a log-normal distribution of width (Te, which is then also the 
width of the Eddington ratio distribution at this fixed true BH 
mass^. And to go from the central to the observed emission- 
line FWHM, we add a scatter to each simulated object se- 
lected from a log-normal distribution of width ctfhwm- Given 
the observed quantities, we can simulate an "observed" virial 
black hole mass. By construction, then, at a given true mass, 
the observed virial mass is unbiased in the mean, and the un- 

^ Here we have assumed that bolometric luminosity is always perfectly 
proportional to continuum luminosity, because this is how we determined the 
obserx'ed bolometric luminosities in previous sections. 
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Fig. 11. — Distribution of quasars from the whole sample in the mass-luminosity diagram, where different colors show BH masses using different virial 
estimators (red for H/3, green for Mgll and blue for CIV). The three diagonal lines (from top to bottom) show 0.01, 0.1 and LEdd respectively. Most quasars are 
bounded by Lbol/^Edd =0.01 and 1. 

in eqn. (|2]l, and the magnitude of the uncorrelated scatters in 
both FWHM (ctfwhm) and Eddington ratio (cte), as well as the 
slope of the underlying power-law BH mass distribution. 

The observed quasar luminosity function has a power-law 
form at the bright end with slope 7^, ^ -2 (Richards et al. 
2006a; Hopkins et al. 2007). We therefore start by assum- 
ing a power-law BH mass function with slope 7^- The slope 
of the simulated bolometric luminosity distribution will be 
^ 7m/(1 + C2) in the limit of small <te, which is required to 
match with 7^, ^ -2. Although there is an observed break in 
the luminosity function (and possibly in the mass function as 
well), we do not include such a break in our modeling, but 
return to that possibility below. 



certainty of virial estimator at this fixed true mass is: 

(Tvir = V {bcF-Ef- + (2ctfWHm)^ ■ (4) 

The two scatters we have added, cte and ctfwhm, are sta- 
tistically independent. We discuss the physical interpretation 
of these scatters, the effects of any possible correlated terms, 
and how these quantities relate to the physical quantities of 
true BLR size and BLR virial velocity at fixed true BH mass, 
in §1531 

Once we have simulated luminosities for each object we can 
impose various luminosity cuts and make comparisons with 
observations. We have the freedom to vary the normalization 
(C\ ) and mass dependence (C2) of the central Eddington ratio 
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Fig. 12. — Apparent distributions of Eddington ratios log(Lbol/i'Edd) based on virial BH masses in various redshift and luminosity bins. Also plotted are the 
fitted Gaussians, with the means and dispersions shown in the middle-left in each panel. The number of objects in each bin is shown in the upper-right comer. 
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TABLE 2 
Model parameters 



Redshift 


7M 


Ci 


C2 




fFWHM 












(dex) 


(dex) 


(dex) 


OJ <z< 1.0 


-2.6 


-3.60 


0.3 


0.4 


0.11 


0.33 


1.9<z<2.1 


-2.6 


-2.88 


0.2 


0.4 


0.11 


0.31 



Note. — See the text for symbol meanings. 

The parameters Ci, C2, and ctfwhm are jointly constrained 
by the observed FWHM distribution in our sample: if the cen- 
tral Eddington ratio is too large or too small, the central value 
of FWHM will be correspondingly under(over)estimated; if 
the slope C2 between logMBH.tme and (logLboi/i^Edd) is too 
small, the central FWHM will show a systematic trend with 
BH mass and luminosity stronger than the one we actually 
observe; if (Tfwhm is too large, then the observed FWHM 
distribution will be too broad. The comparison between the 
observed FWHM distribution and the simulated one must be 
made in each luminosity bin separately. Once we are able to 
reproduce the FWHM distribution and the luminosity distri- 
bution, we automatically reproduce the virial BH mass and 
apparent Eddington ratio distributions. 

Finally, cte and ctfwhm are jointly constrained by the re- 
quirement fjvir = \/ ibcTE)^ + (2(Tfwhm)^ ^ 0.3 dex. If both cte 
and (Tfwhm are too small, we will predict an uncertainty of 
(Tvir < 0.3 dex, an accuracy below that of rever beration map- 
ping. We will discuss this point further in 35.31 

4.4. Results 

We carry out the comparison between our simulation for 
a given set of model parameters and the observed distribu- 
tions at three luminosity bins in each redshift bin. We sys- 
tematically and manually search the parameter space to find 
values where the simulated distributions agree with observa- 
tions, judged by eye. The best-fit parameters are listed in Ta- 
ble |2] We have not searched the parameter space exhaustively 
enough to claim unequivocally that these solutions are unique, 
but the arguments we gave at the end of the previous section 
suggest it is unlikely that there is another region of param- 
eter space which will satisfy all the constraints. Of course, 
more complicated models (e.g., deviations from a pure power 
law in the mass distribution) could also be found with differ- 
ent parameters. Note also that we do not quote formal errors 
on the parameters, given the way in which we found the best 
solution. 

Given the sets of parameters in Table |2] eq. |4] yields an 
uncertainty of the virial masses at fixed true BH mass of 
tJvir = 0.33 dex and 0.31 dex for the Mgii and Civ cases re- 
spectively**. The central Eddington ratio for a typical quasar 
with MBH,true = 10** Mq is (log Lboi/iEdd) = -1-2 and -1.3 for 
the two cases; the (insignificant) difference arises because 
of the difference in the virial relations and in the observed 
FWHM distributions. These distributions of Eddington ra- 
tios must be interpreted with care, since they are heavily de- 
pendent on the exact forms and scatter of the virial relations, 
which are not well-understood at this time. 

The comparison of the simulated and observed results are 
displayed in three luminosity bins in Figs. [T3] (for Mgii) 

* These values are slightly smaller than the observed values (e.g., McLure 
& Jarvis 2002; Vestergaard & Peterson 2006), but our tTvii is the scatter at 
fixed BH mass, while the quoted uncertainties are determined from a sample 
of reverberation mapping objects covering a range of BH masses. 
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Fig. 13. — Comparison between observed and simulated distributions for 
Mgll at 0.7 < z< 1.0. The assumed underlying true BH mass distribution is 
a power-law with slope ■jm = -2.6. The assumed Eddington ratio distribution 
at fixed true BH mass is logLbol/^Edd = ~3.6 + 0.31ogMBH.tiuc + crE, where 
(TE = 0.4 dex. The survey is complete above log = 45.9, and the ob- 

served slope 7l ~ —2 in the power-law distribution of bolometric luminosities 
is reproduced. Each comparison is done for three luminosity bins arranged 
in rows. In each panel, the solid histograms are observed distributions; the 
dotted histograms are simulated distributions for virial masses and apparent 
Eddington ratios based on virial masses. In the BH mass and Eddington ratio 
panels, we also show the distributions of the true underlying BH masses and 
Eddington ratios in dashed histograms. The simulated FWHM distributions 
in the third column are shown as dashed histograms. In the second and third 
column, the numbers show the mean and a of the fitted log-normal for both 
the observed (upper) and simulated (lower) distributions. 
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Fig. 14. — Comparison between observed and simulated distributions for 
CIV at 1.9 < z < 2.1. The assumed underlying true BH mass distribution is a 
power-law with slope 7m = -2.6. The assumed Eddington ratio distribution at 
fixed true BH mass is logLboi/iEdd = ~2.88 +0.2 log Mbh, true + o"E with cte = 
0.4 dex. The survey is complete above log ^ = 46.6, and the observed 

slope 7/. ~ —2 in the power-law distrib utio n of bolometric luminosities is 
reproduced. Notations are similar to Fig. 1131 
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and [T4l (for Civ). The figures show the distributions of BH 
masses (left), Eddington ratios (middle) and FWHMs (right). 
Our simulations are very successful at reproducing the ob- 
served distributions of virial BH masses, Eddington ratios 
and FWHMs. The virial BH masses (solid and dotted his- 
tograms) on average overestimate the true BH masses (dashed 
histograms) in each luminosity bin because of the Malmquist 
bias caused by the finite scatter of the black hole mass esti- 
mator (the part of scatter in virial masses that originates from 
luminosity) and the fact that there are more low-mass than 
high-mass black holes. 

The widths of the observed BH mass and Eddington ratio 
distributions in each luminosity bin are smaller than those of 
the true quantities, and the Eddington ratios are biased low. 
However, although there are more quasars at the faint lumi- 
nosity end, these faint objects are rarely scattered into high- 
luminosity bins, because of the finite width of the true Ed- 
dington ratio distribution. Running our model in alternate 
redshift bins (1.5 < z < 1.7 for Mgii and 3.0 < z < 3.2 for 
Civ) yields similar parameters. Thus our model is consistent 
with a non-evolving true Eddington ratio distribution, at least 
at the current levels of observations. 

There are objects with Lboi/i^Edd > 1 in Figs. [T3]andfT4l 
caused by the high-luminosity tail in the scatter cte = 0.4 dex at 
individual BH masses. However, the fraction of these super- 
Eddington objects is small and the formal Eddington limit is 
only approximate. Enforcing a cutoff in the model Eddington 
ratio distribution has little effect on our conclusions. 

Malmquist bias will be present independent of any luminos- 
ity cut on the sample, given the scatter between the true and 
"observed" masses and a bottom-heavy true BH mass distri- 
bution. This Malmquist-type bias is important in analyses that 
are based on "observed" BH masses, such as the active BH 
mass function (e.g., McLure & Dunlop 2004) and the redshift- 
evolution of the M-a relation (Lauer et al. 2007). 

We have demonstrated the Malmquist bias with simula- 
tions, but the form of the bias can be shown analytically. We 
start by deriving the bias at fixed observed virial mass, and 
then generalize to the bias at fixed observed luminosity. Fol- 
lowing Lynden-Bell et al. (1988), we model the probability 
distribution of m^, = InMBH.vh- given m = InMeH.true as a Gaus- 
sian; that is: 



PQ{me\m) = {liT Al^) '/^exp 



(5) 



where A,„ = (InlO)crvir is the Gaussian a. Let niilnMBH,tiue 
be the number of BHs with masses between InMBH.tme and 
InMBH.true + t^lnAfBH.tme- For a power-law distribution of the 
underlying true BH masses, n oc M2 



^BRtrue'Wehave: 



dMu 



ndm^ MZ L dMmm. « e"""" dm , (6) 



^BH.tme 

where jm is our assumed slope in the underlying power-law 
BH mass distribution. Then the distribution of m for a given 
is 

^ -1 

pi(m\m,) = poe''"'" I poe^"'"dm ' 



: (27r )-i/2 exp |-i [m - (m, + 7m A^ )] V A^ 



(7) 



This is a Gaussian distribution, centered on the mass 
(m) = me + 7m A,^, and therefore the bias in log mass is 



-7MA^,/lnlO. Inserting 7m = -2.6, cTyir ~ 0.3 and A,„ = 
(In 10)(Tvir ~ 0.7 produces a bias in log mass of 0.55 dex. 

However, the Malmquist bias we derived above is the bias 
at fixed virial mass, not the bias at fixed luminosity (Figs. [T3] 
andfT4li. To calculate the latter, we proceed similarly. 

Let / = InLboi and m = InMBH.tme, then the probability dis- 
tribution of / at fixed m is: 



poa|m) = (27rA^)-'/2exp 



[l-iC-i + CAm)] 



21 



2Aj 



(8) 



where the constants C3 = [Ci +log(l .3 x 10^^)] In 10, C4 = 1 + 
C2 and A/ = (In 10)(7e following equation (|2]i. 
The probability distribution of m at fixed I is then 



pi(m\l) = Poe''"'" 



poe'"^"'dm 



-1/2 



exp 



'-4 

2Aj 



CaH- 



-C3) + 7mAP 



<-4 



(9) 

which is a Gaussian distribution with dispersion A/ /C4 cen- 
tered on (m) = (/ - C3)/C4 + 7mA^/C|. A similar analy- 
sis for the distribution of the "observed" BH mass = 
InMflH.vir at fixed luminosity gives a Gaussian distribution 
with mean (m^.) = (/-C3)/C4 and dispersion 2crFWHM InlO. 
The Malmquist bias, i.e., the difference in the two means is 
-7mA^ / (C| In 10) « 0.6 dex when plugging in numbers. The 
distribution of log true masses at fixed luminosity has a width 
of A/ / (C4 In 10) = 0.32 dex; while the distribution of log virial 
masses at fixed luminosity has a width of 0-^^. = 2 (Tfwhm = 
0.22 dex. These results are in excellent agreement with what 
we see in Figs. [T3]and[T4l 

The Malmquist bias is mainly caused by the bottom-heavy 
true BH mass distribution and the scatter in the virial es- 
timators (or the uncorrected scatter cte in the case of fixed lu- 
minosity). However, we emphasize that the exact magnitude 
and form of the bias depends in detail on the partition and 
behavior of the scatter in the virial relations. We have used 
(Je = 0.4 dex and ctfwhm = 0.11 dex in order to reproduce the 
observed uncertainty in virial BH estimators of dvir ?J 0.3 dex 
and the observed FWHM distributions, but a smaller scatter of 
virial BH mass at fixed true mass would reduce the Malmquist 
bias. Also, the values of cte and ctfwhm could depend on lu- 
minosity. Our data do not allow us to consider such refine- 
ments in detail. Since the Malmquist bias is proportional to 
7m, we expect a smaller bias for objects fainter than the break 
in the luminosity function, such as those in the AGES sample 
(Kollmeier et al. 2006). 

Although our model of the underlying true BH mass and 
Eddington ratio distributions can roughly reproduce the ob- 
servations, we have made a few simplifications. In partic- 
ular, the largest uncertainty arises from the assumed virial 
relations and their scatters, for which there is no consensus 
yet. Changes in these virial relations (e.g., a different slope 
in the R-L relation, etc., see Bentz et al. 2006; Netzer et al. 
2007) will certainly change the values of our model parame- 
ters. We also note that there might be other systematic errors 
associated with virial BH masses or reverberation mapping 
masses in general (Krolik 2001). Given all the limitations and 
simplifications of our model, we do not wish to overinterpret 
our modeled distributions at this time, but rather highlight the 
danger that the virial masses for a population of quasars could 
be subject to serious systematic biases. Until we have both 



16 



SHEN ET AL. 



TABLE 3 

QUASAR BLACK HOLE MASS FUNCTION 



z 




*Mbh>io« Mq (Mpc-3) 


0.85 


2.87 X 10"'' 


1.0 X 10-** 


2.0 


1.20 X 10-5 


7.4 X lO-'* 


3.1 


8.08 X 10"* 


5.0 X lO-** 



better observations and a theoretical understanding of the sys- 
tematic behavior of the virial estimators with BH properties, 
we urge caution in the interpretation of BH mass functions 
based on virial estimators. 



4.5. Active BH mass function in quasars 

Our model explicitly provides the true BH mass distribution 
and connects it to the bolometric luminosity distribution via 
our modelled Eddington ratio distribution, and thus provides 
an estimate of the underlying BH mass function in quasars 
based on our knowledge of the bolometric luminosity func- 
tion. 

An alternative method using optically selected samples uses 
the observed luminosity function and the fraction of BHs 
above a certain mass threshold within some luminosity range, 
to place lower limits on the number density of active BHs 
above the mass threshold (McLure & Dunlop 2004). The 
measurement is a lower limit simply because it neglects BHs 
scattered out of the luminosity range being considered. Using 
this method, McLure & Dunlop (2004) found that the majority 
of > 10*^^ Mq SMBHs are already in place at z 2 when geo- 
metric obscuration and quasar duty cycles were taken into ac- 
count. However, their results were based on virial BH masses, 
and they did not take into account the Malmquist bias that we 
have discussed here. 

Here we use our model to estimate the cumulative BH mass 
function in quasars above a certain mass threshold. Notice 
that the break in the bolometric luminosity function (Hopkins 

et al. 2007) is at log (^^^pr) 46 at z 1 - 3, correspond- 
ing to logMBH.true ^ 8.8-9.0 in our model (e.g., equation|2]l. 
Thus our single power-law distribution of the underlying BH 
masses is inaccurate below this mass threshold^, and will af- 
fect the mapping from the BH mass function to the bolometric 
luminosity function. We will explore the effects of a break in 
a future paper. 

We focus on the cumulative mass function above 10^'^ Mq, 
i.e., the most massive SMBHs. The estimation is simple: we 
take our simulated catalogs at different redshifts and count 
the number of BHs with mass greater than 10^'^ Mq and 

the number of BHs with log (^^^j^^ > 46; the latter number 

is compared to the cumulative bolometric luminosity func- 
tion (Hopkins et al. 2007, the "full" model in their table 3) 
with logLboi > 46 and normalized; then we use this normal- 
ization to estimate the number density of active BHs with 
logMBHtme > 9.5. Our results are listed in Table[3] for three 
redshift bins: z = 0.7- 1.0, 1.9-2.1, 3.0-3.2. The number 
density of BHs with mass > 10*^^ Mq in optical quasars peaks 
around z^2, and then decreases at lower redshift. 

The numbers above are lower limits on the total (active 
and inactive) SMBH population because of quasar duty cy- 
cles and the presence of obscured populations of active BHs 

' Our conclusions on the distributions of FWHMs and virial BH masses 
are insensitive to this detail. However, the Malmquist bias will be smaller 
below the break, as mentioned above. 



at high redshift. Surveys in other wave bands, especially in 
the X-ray and infrared, are starting to reveal those missing ac- 
tive SMBHs (e.g., Treister et al. 2006; see also Zakamska et 
al. 2003), and will provide better constraints of the BH mass 
function at high redshift, given our knowledge of the local 
dormant BH mass function (e.g., Marconi et al. 2004; Mer- 
loni 2004). We note that the numbers derived here are fully 
consistent with current constraints on the local BH mass func- 
tion at the high-mass end (e.g., Tundo et al. 2007; Lauer et al. 
2007), and allow an additional boost factor of ^ 10-20 due 
to quasar duty cycles and geometric obscuration at z ^ 2. 

4.6. Results for radio loud quasars and BALs 

In this section, we determine virial BH masses for radio- 
loud quasars and BALs, to explore whether their black hole 
properties are distinctive from "ordinary" quasars of similar 
luminosity and redshift. The Eddington ratio and BH mass 
distributions will be affected by Malmquist bias, but the bias 
should be similar for each individual subsample, allowing this 
kind of comparison to be made. 

For the radio analysis, we carve out a narrow range in 

redshift-luminosity space, log (^j^^r^ = 46.6-46.75 and z = 

1.5-2.3, over which we have Mgii virial masses for each 
source, and the FIRST-detected and undetected quasars have 
similar redshift and optical luminosity distributions. At these 
redshifts, the FIRST radio flux limit corresponds to a radio 
luminosity density of 3 x 10^^ Watts/Hz, which is radio- 
loud by all standard definitions. This region contains 3360 
radio-undetected quasars and 170 radio-loud quasars. The 
radio-loudness, defined as TZ = f(, cm/fim) for these FIRST- 
detected objects ranges from ~ 5 to ^ 4600 with median value 
^ 100, based on a subset of DR3 quasars studied by Jiang et 
al. (2007b). Fig. [15] shows that the Mgii based virial mass 
distributions for the radio-undetected and radio-loud quasars 
are slightly different, with radio-loud quasars having slightly 
larger median value by ^ 0.12 dex, comparable to the results 
in McLure & Jarvis (2004). Radio-loud quasars in our sample 
have smaller Eddington ratios than radio-undetected quasars 
on average, as has been found in previous analyses (e.g.. Ho 
2002). 

For BALs (most of which are Civ BALs), the situation is 
more complicated. First, BALs in general are more reddened 
than nonBALs (e.g., Reichard et al. 2003), affecting the bolo- 
metric luminosity estimates. Second, the broad absorption 
troughs can severely impact the reliability of FWHM mea- 
surements; the concerns about disk winds affecting the mea- 
surement of the Civ line (discussed in § 13.3.21 ) are especially 
important for BALs. 

Fortunately, for 1.7 < z < 2.2 we have both Mgii and Civ, 
thus it becomes feasible to compare the difference in virial 
mass for high-ionization BALs (which appear in the Civ line 
only) and nonBALs, based on the Mgii estimator For the 
uniformly-selected quasar sample the redshift and luminosity 
distributions are almost identical for Civ BALs and nonBALs 
at 1-7 < z < 2.2 (Shen et al. 2008). We choose uniformly- 
selected quasars that have redshift 1.7 < z < 2.2 and have 
Mgii based virial masses; the sample includes 5408 nonBALs 
and 796 BALs. Fig.[T6lshows that their Mgii virial mass dis- 
tributions are quite similar, with median mass logMBH.vii- = 
9.34 for both. This result is consistent with recent findings 
by Ganguly et al. (2007) based on 5088 quasars from SDSS 
DR2, although the mean virial BH masses are different due 
to difference in the calibrations used and line measurements. 
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Fig. 15. — Distributions of Mgll based virial masses for tlie radio- 
undetected (solid histogram) and radio-loud (dashed histogram) quasars in 
subsamples matched in redshift and luminosity distributions; see text. 
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Fig. 16. — Distributions of Mgll based virial masses for nonBALs (solid 
histogram) and BALs (dashed histogram) with 1.7 < ; < 2.2, which are 
within the uniformly-selected quasar sample. The redshift and luminosity 
distributions of nonBALs and BALs are identical. 



as we discuss below. This similarity in virial BH masses for 
BALs and nonBALs supports the idea that BALs are intrinsi- 
cally no different from nonBALs, and the BAL phenomenon 
is caused by different viewing angles with respect to a disk 
wind (e.g., Weymann et al. 1991; Elvis 2000). This picture is 
also consistent with our recent findings that BALs show sim- 
ilar large-scale clustering to nonBALs (Shen et al. 2008). 

5. DISCUSSION 

5.1. General issues with virial estimators 

The virial BH mass measurements presented in this paper 
are consistent with, but not identical to results from a variety 
of papers (e.g., McLure & Dunlop 2004; Dietrich & Hamann 



2004; Kollmeier et al. 2006; Ganguly et al. 2007). The differ- 
ence in the mean of the distribution can be due to a variety of 
effects: different samples, different virial estimators and dif- 
ferent calibrations for the same estimator, as well as different 
procedures to measure luminosities and line widths used in 
these virial estimators. We elaborate below how these various 
choices lead to differences in virial masses even for a given 
underlying sample. 

5.1.1. Different virial calibrations for the same line estimator. 

There are several versions in the literature of a given virial 
estimator, and they do not generally yield the same virial mass 
when applied to SDSS or other quasar/ AGN samples. For ex- 
ample, the H/3 calibrations in McLure & Dunlop (2004) and 
in Vestergaard & Peterson (2006) differ by ~ 0.17 dex given 
the mean luminosity ALa,5ioo ~ lO^"^'^ erg s~' in our sample. 
Likewise, for the Mgii estimator, the version in McLure & 
Dunlop (2004) yields a larger virial mass by ^ 0.3 dex for 
AiA.3000 = 10"*^ ' erg s"' (which corresponds to SDSS quasars 
at z ~ 2), than the older version in McLure & Jarvis (2002). 
The latter case partially explains the discrepancy between our 
measurements and a subset of our sample measured by Gan- 
guly et al. (2007). While we have used the most recent cal- 
ibrations available (McLure & Dunlop 2004; Vestergaard & 
Peterson 2006), it is quite likely that those calibrations will 
change when the updated R-L relation and virial coefficient 
(zero-point offset) are incorporated (Onken et al. 2004; Bentz 
et al. 2006; Kaspi et al. 2007). 

5.1.2. Different ways to measure luminosity and FWHM. 

Even for the same calibration, different authors have used a 
variety of ways to measure luminosity and line width param- 
eters, which will lead to systematically different results. For 
example, the mean Mgii FWHM taken from the broad com- 
ponent of a two-Gaussian fit is larger by 0.15 dex than the 
non-parametric FWHM measured from the modelled spec- 
trum which best reproduces the line profile. This contributes 
to the underestimation of the Mgii based virial masses (e.g., 
Dietrich & Hamann 2004; Ganguly et al. 2007). For Civ, 
on the other hand, it is more reasonable to use the FWHM 
measured from the modelled flux that best reproduces the Civ 
line profile because that is how the Civ relation is calibrated 
(e.g., Vestergaard & Peterson 2006). Although some argue 
that a narrow-line component is present in the Civ line (e.g., 
Marziani et al. 1996) and recipes have been provided for sub- 
tracting it (cf., Bachev et al. 2004), we cannot perform such an 
exercise until the Civ estimator is re-calibrated with narrow- 
line-subtracted line widths. Indeed, it is crucial to use the 
methods of measuring FWHMs used in the original calibra- 
tions of these virial relations if one wants to obtain unbiased 
results. We feel confident about our choices, as we found 
no mean offset between virial masses based on different lines 
(Figs.|5]and|9]l. 

5.1.3. Other issues with our measurements. 

There are some additional factors that might lead to dif- 
ferences between our measurement and others for the same 
sample of quasars. First, the spectroscopic flux scale used in 
this study is higher than previous SDSS data releases (DR5 
and prior) by 0.14 dex, giving rise to a systematic ^ 0.07 
dex increase in virial masses. Second, we have separated the 
broad and narrow components of H/3 and Mgii at a FWHM of 
1200 km s"', following Hao et al. (2005). It is possible that 
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we are biasing the mean virial mass by doing so. Neverthe- 
less, these effects are minor considering the uncertainties in 
the zero-point offsets of these virial estimators. 

As we discussed in 34.21 given the independence of the 
observed FWHM distribution with luminosity, objects in the 
A^BH.vir - Lbo\ plane tend to lie on a narrow stripe with a 
slope b ^ 0.5 (see Fig. [TTT i. The lower-luminosity high red- 
shift quasars in the Kollmeier et al. (2006) sample would 
also fall on the stripe. The mean apparent Eddington ra- 
tio (log(Lboi/iEdd))vii- oc (1 -/7)logLboi- Thus the expected 
(log(Lboi/LEdd))vir would be ~ -1.75 for log (^) =44, ^ 

-2.25 for log (^) = 43, and ^ -3.25 for log (^) = 41 

for b ~ 0.5. These Eddington ratio values are in good agree- 
ment with the typical Eddington ratios from X-ray selected 
AGN samples which cover the bolometric luminosity range 
iboi ~ 10**' - 10'^ erg s-i (e.g., Panessa et al. 2006; Babic et 
al. 2007; Ballo et al. 2007). However, there are two reasons 
to be cautious about this apparent agreement: current X-ray 
selected AGN samples still suffer from limited sample size 
and various incompletenesses, and most of the BH masses 
in X-ray selected AGN are estimated based on host galaxy 
properties rather than virial methods, and it remains unclear 
how good these methods agree with each other for those low- 
luminosity BHs. 

On the other hand, virial methods have their own limitations 
due to their reliance on poorly understood BLR physics. If the 
existence of the broad line region itself depends on Eddington 
ratio, then it is not surprising that broad line AGN/quasars lie 
preferentially within a favorable Eddington ratio range. More- 
over, the width of the Eddington ratio distribution at fixed lu- 
minosity could be underestimated, if, for example, the line 
widths only partially reflect the v irial velocities in the BLR, 
as we wiU further discuss in 35.31 

5.2. Limitations of the SDSS sample 

Although the large size of the SDSS sample provides un- 
precedented statistics, we feel obligated to point out its limi- 
tations. 

First, the dynamic range in luminosity is narrow for SDSS 
quasars. Only luminous quasars are included in the flux- 
limited sample. Other surveys, such as the AGES survey 
(Kochanek et al. 2004) and the 2dF-SDSS LRG and QSO 
(2SLAQ) Survey (Cannon et al. 2006) extend to fainter lumi- 
nosities at the same redshift. Our sample does include some 
objects of lower luminosity, but they are incomplete in ways 
that are hard to quantify. 

Second, the wavelength range of SDSS spectra has forced 
us to use different virial estimators at different redshift. Al- 
though we cross-compared virial relations of H/3-Mgii and 
Mgii-Civ, the direct comparison between the H/3 and Civ 
estimators is impossible with SDSS spectra. Thus it is im- 
portant to study these broad lines simultaneously with multi- 
band spectra (e.g., Sulentic et al. 2000 and references therein; 
Baskin & Laor 2005). It is also desirable to obtain near-lR 
spectra for high-redshift quasars with optical spectra, to get 
H/3 or Mgii measurements (e.g., Sulentic et al. 2006; Jiang et 
al. 2007a; Kurk et al. 2007; Netzer et al. 2007), since the Civ 
virial estimator alone is perhaps questionable for the reasons 
we have discussed in this paper. 

5.3. Issues with our model 



In 34.31 we introduced a statistical model which can repro- 
duce the observed luminosity and FWHM distributions. We 
now provide some justifications of our choices of model pa- 
rameters, and discuss the connections between the real physi- 
cal quantities, i.e., BLR size and virial velocity, and their sur- 
rogates, luminosity and FWHM. 

An important assumption of our model is the origin of 
the uncertainties in virial estimators, i.e., for fixed true BH 
mass, the virial estimator will give an estimate log-normally 
distributed around the right mean value and with dispersion 
(Tvii. However, the value of cTvii- remains unclear: the rms 
scatter of virial masses around reverberation mapping (RM) 
masses (e.g., McLure & Jarvis 2002; Vestergaard & Peter- 
son 2006) is 0.3-0.4 dex, comparable to the rms scatter of 
virial masses around BH masses derived from the M-a re- 
lation (e.g., Greene & Ho 2006); but the relationship of this 
number to the scatter m'ound fixed true BH mass is unknown. 
Given all the issues with virial estimators we discussed in 
35.11 it seems unlikely that virial estimates based on single- 
epoch spectra are good to a factor of 2, an accuracy that can 
barely be achieved with reverberation mapping. Therefore we 
believe our choice of cTvii ^ 0.3 dex is appropriate. 

At a given true BH mass, in our simplistic model, the broad- 
ening of the luminosity distribution is completely uncorre- 
cted with the broadening of the FWHM distribution. Are our 
required values for the uncorrelated scatters cte and ctfwhm 
consistent with observations? Current reverberation mapping 
data indicate a scatter of ~ 0.2-0.35 dex in luminosity at 
fixed BLR size (e.g., Kaspi et al. 2005; Bentz et al. 2006). 
The best studied reverberation mapping sources indicate scat- 
ters of ^0.1 dex in line width at fixed BLR size for a given 
object (e.g., Bentz et al. 2007). These values are for a handful 
of well studied objects with reverberation mapping of the H/3 
line only, and we might expect somewhat larger scatters for 
both luminosity and line width in samples with a wider range 
of luminosity and redshift than probed by current RM sam- 
ples. Thus while our choices of uncorrelated scatters cte = 0.4 
dex and (Tfwhm = 0.1 1 dex may appear rather large, they are 
at least plausible. 

We still expect some component of the variations in lumi- 
nosity to correlate directly with variations in FWHM, as is 
seen in local samples. Here we investigate what magnitude 
of correlated broadening can be supported given existing ob- 
servations. At fixed true BH mass, the virial mass estimate is 
expressed in terms of its luminosity and FWHM: 

logMfiavn = C+felogLboi + 21ogFWHM 

= C + fe((logLboi) + (5E + 4o,T) (10) 

+ 2((logFWHM) + 5FWHM-0.5Mco,T) , 

where b is the slope of the R-L relation, the various (5,- 
denote Gaussian-random variables with dispersions ct,-, and 
other constants have been absorbed in C. At this fixed true BH 
mass, the bolometric luminosity and the FWHM both follow 
a log-normal distribution, logLboi = (logLboi) + (^e + (5con- and 
log FWHM = (log FWHM) + (5fwhm- 0.5 WcoiT- In addition to 
the uncorrelated terms 5^ and (5fwhm as modelled in 34.31 we 
also introduced correlated variations 5^o„ and -Q.5b5co„ in lu- 
minosity and FWHM, respectively, whose amplitude we will 
constrain below. Note that these correlated terms cancel each 
other out, and hence do not contribute to the uncertainty in the 
virial relation. 

The correlated terms in equation (fTOb reflect the intrinsic 
variation in BLR size/virial velocity at fixed true BH mass. 
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Fig. 17. — Schematic plot showing the relations between the actual phys- 
ical quantities, i.e., BLR radius logR and virial velocity logV that deter- 
mine the BH mass, and their surrogates, luminosity logL and line width 
logFWHM. The mean R — L relation is assumed to have a slope 0.5. The 
mean relation between logFWHM and log/? is assumed to have a slope —0.3 
in this particular example, i.e., not a perfect virial relation with slope -0.5. 
The intrinsic variations of BLR radius and virial velocity are 5 and -0.55 at 
this fixed true BH mass. For luminosity, 5i is the rms scatter around the R—L 
relation and 82 = 25 is the variation that drives the variation in BLR size. For 
FWHM, 5fwhm is the rms scatter around the mean FWHM-virial velocity 
relation, and ^3 = -0.35 is the part that responses to the variations in BLR 
size. Terms in red end up to be the uncorrelated variations (5e and <5fwhm in 
our model, which contribute to the uncertainty of the virial estimator Terms 
in blue end up to be the correlated variations that do not contribute to the 
virial uncertainty. 



In our model formalism in §43] we have set CTcoit = 0. Since 
our observed FWHM distributions are quite narrow, we found 
that we need'" (icon- ^ 0.2 dex, such that the distribution of 
simulated FWHMs will not be broader than the observed dis- 
tribution, keeping the requirement that Uvir > 0.3 dex. With 
this limit, the intrinsic dispersion in BLR size for our sam- 
ple is (Tblr 01 dex at fixed true BH mass, which is nar- 
rower than the typical dispersion in BLR size inferred from 
less luminous reverberation mapping sources. It is possible 
that the intrinsic dispersion of BLR size at fixed true BH mass 
is indeed smaller for higher-luminosity objects. Reverbera- 
tion mapping data for high luminosity objects are required to 
test this scenario. 

One could also get a small apparent dispersion in BLR size 
if systematics in the BLR virial velocity-line width relation 
cause the FWHM to only partially trace real variations in 
the virial velocity. In other words, the mean relation in the 
log FWHM -log/? diagram for fixed BH mass could have a 
slope shallower than -0.5, as has been suggested by Bentz et 
al. (2007) in H/3 reverberation studies in NGC 5548. A non- 
virial-gas contaminated line width and/or inappropriate meth- 
ods to measure line width could lead to this problem. The 
effects of this are illustrated in a specific worked example in 
Fig. [17] which shows how the intrinsic variations in the true 
BLR size {&) and in the true virial velocity (-0.5(5, exactly 
compensating) for a given true BH mass might be parcelled 
out. The dispersion of observed luminosity is composed of 

This upper limit is set by the condition that the width of the correlated 
part -0.5/)5coiT is half of ctfwhm in FWHM so that its contribution to the 
broadening of the FWHM distribution is negligible when added in quadra- 
ture. 



the portion that fully reflects the intrinsic variation in BLR 
size 62 = 2S (where b = 0.5 is assumed), and the portion of 
rms scatter around the mean R-L relation Si. The disper- 
sion of FWHM is composed of the portion that only partly 
reflects the intrinsic variation in virial velocity -0.3d (i.e., a 
slope of -0.3 is assumed for the mean log/?-logFWHM rela- 
tion in this particular example), and the portion of rms scatter 
around the mean virial velocity-FWHM relation, (5fwhm- Di- 
viding 62 into two terms: the Scon = 1-2(5 term cancels with 
the -0.3(5 variation in FWHM, i.e., these are the correlated 
terms defined earlier; the remaining 0.8(5 portion in 62 com- 
bines with (5i to form 5e. Hence, in this case, although Se is 
the uncorrelated term in luminosity according to our original 
definition, it still includes a portion that reflects the intrinsic 
variation in BLR size, and it contributes to the uncertainty in 
virial estimators. The intrinsic variation in BLR size 5 can 
now be slightly broader, i.e., (Tblr ^ 0.17 dex. 

6. CONCLUSIONS 

We have measured virial BH masses for 58,664 quasars in 
the SDSS DR5 quasar catalog. We used and compared three 
virial estimators: H/3 (z < 0.7), Mgii (0.7 < z < 1.9) and Civ 
(z > 1 .9). We emphasized the importance of using the original 
definitions of line width and luminosity in whichever virial 
calibration is used. Our main conclusions are the following: 

1 . Within our sample, the line widths follow a log-normal 
distribution; their means and dispersions depend only 
weakly on redshift and luminosity. 

2. For a subsample of quasars for which both the H/3 
and the Mgii estimators are available, the ratio of their 
FWHMs follows a log-normal distribution with mean 
0.0062 and a dispersion of 0.11 dex; the ratio of virial 
BH masses based on the two lines also follows a log- 
normal with mean 0.034 and a dispersion of 0.22 dex. 
Therefore, the Mgii and H/3 estimators give consistent 
results. 

3. We further compared the Mgii and Civ estimators in a 
subsample of quasars with both lines. Their FWHM 
ratio follows a log-normal with mean 0.027 dex and 
dispersion 0.18 dex; the ratio of their virial estimates 
follows a log-normal with mean -0.06 dex and dis- 
persion 0.34 dex. Thus virial BH estimates are con- 
sistent using both lines, although with larger scatter 
than between H/3 and Mgii. However, the dispersion 
in the FWHM ratios is comparable to or even larger 
than the dispersion of the Mgii and Civ FWHM dis- 
tributions themselves, and the Mgii FWHM and Civ 
FWHM are weakly correlated at best (see Fig.[8]l. Fur- 
thermore, the Civ estimator tends to give smaller virial 
masses than the Mgii estimator for objects with small 
blueshifts (< 1000 km s"'), and conversely for objects 
with larger blueshifts. The reason for this systematic 
bias is attributed to the geometry of the BLR. The Civ 
line shows many features that suggest a component 
from non-virialized gas such as a disk wind (Murray 
et al. 1995; Proga et al. 2000; Elvis 2000). Therefore 
the Civ FWHM is perhaps not a good indicator of the 
BLR virial velocity. The current calibration of the Civ 
estimator gives consistent results with those using the 
other two estimators in the mean, but we caution that 
the bias may be large for individual objects. 
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4. The typical range of virial BH masses in the SDSS 
quasar sample is 10^- 10'° Mq. The upper envelope 
of the virial mass distribution rises up to z 2 and then 
flattens out. There is a clear upper limit ^ 10'" M© for 
all quasars, as other studies have found (e.g., McLure 
& Dunlop 2004; Vestergaard 2004). 

5. Quasars lie in a narrow stripe in the mass and luminos- 
ity diagram (Fig. fTTT i bounded by Lboi = O.OlLEdd and 
Lboi = i'Edd- consistent with recent findings (e.g.. Woo 
& Urry 2002; Kollmeier et al. 2006). However, this 
distribution is implicitly constrained by the virial rela- 
tions and the observed FWHM distributions. Similar 
to the findings of Kollmeier et al. (2006), the distribu- 
tions of apparent Eddington ratios Lboi/^Edd based on 
virial BH masses in different redshift-luminosity bins 
follow log-normal distributions, with means between 
Lho\/LEdd « 10"" and 10"°^. The widths of these dis- 
tributions are typically ~ 0.3 dex or less. 

6. The narrowness in the observed Eddington ratio (and 
virial mass) distributions within each luminosity bin is 
interpreted as arising from the combination of luminos- 
ity cuts and the underlying distributions of FWHMs; it 
is not the same as the intrinsic uncertainty in the virial 
mass estimators. By assuming an underlying true BH 
mass distribution and an Eddington ratio distribution at 
fixed true BH mass, together with the assumptions that 
the observed luminosity and FWHM are imperfect trac- 
ers of the virial BH mass, we were able to reproduce the 
observed distributions of luminosities, FWHMs, virial 
BH masses and apparent Eddington ratios in each lu- 
minosity bin. Monte Carlo simulations demonstrate 
that the observed virial BH mass and apparent Edding- 
ton ratio distributions are subject to Malmquist bias, 
i.e., more lower-mass BHs are scattered upwards due to 
the scatter between virial mass and true BH mass than 
higher-mass BHs are scattered down. To better quantify 
this Malmquist bias we need a better understanding of 
the form and scatter in the virial relations. 

7. We also compared the distributions of virial BH masses 
for radio-loud quasars and BALs with that of "ordi- 
nary" quasars matched in redshift and luminosity. The 
mean virial mass of radio quasars is ~ 0.12 dex larger 
than that of ordinary quasars, but the mass distribu- 
tion of BALs is indistinguishable from that of ordinary 
quasars. 

With ever larger quasar samples, it has now become fea- 
sible to measure the clustering properties of quasars, and to 
directly test galaxy formation scenarios within the hierarchal 
structure formation framework (e.g., Shen et al. 2007 and ref- 
erence therein). The clustering measured in current quasar 
samples shows only a weak luminosity dependence at z < 2.5 
(e.g., da Angela et al. 2007). Since BH mass is tightly corre- 
lated with bulge properties (e.g., Tremaine et al. 2002), we ex- 
pect a correlation between BH mass and host dark matter halo 
mass (e.g., Ferrarese 2002). The fact that luminosity does not 
strongly correlate with clustering strength seems to indicate 
that the instantaneous quasar luminosity is not a good indica- 
tor of BH mass (Lidz et al. 2006). In fact, both in our model 
and in the data itself we see a significant range in luminosi- 
ties at a fixed BH mass (see Figs. [13] and [l4] l. We plan to 



study quasar clustering as a function of virial BH mass, tak- 
ing care to incorporate the effects of the uncertainties in the 
virial estimators and the Malmquist bias. 

Despite the biases we have identified in this paper, the virial 
estimators are irreplaceable tools for estimating BH masses 
in AGN/quasars. One must simply be very careful when in- 
terpreting these virial masses for individual objects and for 
statistical samples. We point out that our methodology here 
can be applied to future data sets which push to lower lumi- 
nosity AGN/quasars at all redshifts. It will be particularly 
interesting to probe below the break in the quasar luminos- 
ity function, where the Malmquist bias should be smaller and 
where it has been suggested that the Eddington ratio should 
be smaller as well (e.g., Hopkins et al. 2006). This will allow 
us to explore both the low-mass end of the BH mass function, 
and to study the nature of objects with low Eddington ratio at 
high redshifts. In addition, there is a strong need for better 
understanding of the forms and scatter in the virial relations. 
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